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In dam works, air precooling of aggregate is a common and effective method to avoid temperature cracks in con- 


crete structure. In order to offer a reliable design theory for the air precooling process to avoid unreasonable 


energy consumption, the transient heat transfer characteristics of the aggregate are intensively analyzed. The 


combined structure of the aggregate and the interstitial space in the hopper is treated as a porous structure, and the 


space-average method is used to simulate the transient heat transfer process. Simulation results show that size of 


the hopper and the average air velocity in the cross section have great influence on the transient heat transfer 


process of the aggregate, while the porosity in the range of 0.4-0.5 has little influence. Nomograms are abstracted 


from simulation results, and then correlations of the compared excess temperature are precisely fitted to predict 


the air precooling transient heat transfer process of the aggregate. 
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Introduction 


The safety of the dam construction prominently de- 
pends on the temperature control of the concrete structure 
[1]; precooling of aggregate is an effective method to 
avoid temperature cracks, and is widely used in actual 
concrete works. By now, the aggregate precooling pro- 
cess is very critical for safety and economic performance 
of the hydraulic and hydro-power engineering, and the 
research on it is extremely important. 

During the air precooling process, the aggregate 
moves down from the feeding opening to the discharge 
opening in the aggregate hopper, while the cool air flows 
up from the air inlet to the air outlet. The whole counter 
flow process makes the aggregate and the air fully 
exchange heat in the aggregate hopper. In actual concrete 
engineering, analysis of the air-cooled transient heat 
transfer process is based on simplified equivalent con- 
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version of the transient heat transfer characteristics for 
water-cooled concrete [2-4], and to avoid the engineering 
disqualification, excessive design margin is usually 
adopted and then makes high energy consumption. In 
order to offer a reliable design theory and then rea- 
sonably reduce energy consumption, it is urgent to make 
an intensive study of the air-cooled transient heat transfer 
characteristics of the aggregate. 

In this paper, the aggregate and the interstitial space 
are combined as a porous medium, and the transient heat 
transfer process in the aggregate hopper is simulated by 
means of the porous medium model. The simulations 
almost cover all actual working conditions. According to 
the simulation data, the influences of many factors, such 
as the hopper construction, the porosity, the air velocity, 
the cooling time, etc., are analyzed, and the transient 
cooling regularity of the aggregate in the hopper is 
deduced. 
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Nomenclature 
a heat diffusivity, m7/s T aggregate temperature, K 
A fitting coefficients in Eq. (23) Ta air temperature, K 
Ah cross section area of the hopper, m? To initial temperature of the aggregate, K 
B fitting coefficients in Eq. (23) U perimeter of the cross section, m 
Bi Bice aunahes 7 average air flow velocity in cross section, 
m/s 
c specific heat, J/(kg: K) Ua inlet flow velocity, m/s 
C fitting coefficients in Eq. (23) Ux air flow velocity in interstitial space, m/s 
d stone diameter, m Greek symbols 
D diameter of the hopper. m , Sonia heat transfer coefficient, 
W/(m*: K) 
D’ effective diffusion tensor in Eq. (2) E Porosity 
Fo Fourier number A thermal conductivity, W/(m: K) 
Fo* fitting coefficients in Eq. (23) u kinematic viscosity, kg/(m-:s) 
H height of the arrangement unit, Fig. 7 p density (kg /m’) 
kef effective thermal conductivity in Eq. (2) T time, s 
L height between two air flues, m Subscripts 
Pr Prandtl number a air 
QO; convection heat transfer in interstitial space, Fig. 7 e equivalent values 
Q> heat conduction between two adjacent stones, Fig.7 kx interstitial space 
r radius of the touch area, Fig. 7 S stone 
stone radius, m v virtual values 


Re Reynolds number 


Similarity Experiment and Simulation 


The rough spherical packed bed is usually treated as a 
porous medium composed by the rough spherical materi- 
al and the interstitial space in industry [5, 6]. The struc- 
ture composed of the aggregate and the interstitial space 
in the aggregate hopper is similar with that of the packed 
bed, and also can be treated as a porous structure. The 
whole transient heat transfer process of the porous me- 
dium can be simulated by CFD software, which is suita- 
ble for large stacks of particles, and the simulation results 
show good prediction [7, 8]. 


Similarity experiment 

Owing to the testing restriction on the construction 
site, it is difficult to obtain comprehensive and precise 
testing data, so it is necessary to establish a small-scale 
similarity experimental system to availably analyze the 
transient heat transfer characteristics of the aggregate 
during the air cooling process. 

According to the similarity criterion of the heat trans- 
fer model, the Nusselt number is the same when the heat 
transfer phenomenon is similar [9]. The definition of the 
Nusselt number is: Nu=ad/⁄AÀ, where a is the convective 
heat transfer coefficient imposed on the stone surface, 


and mainly depends on the Reynolds number, so with the 
same Re, as well as the same air flow velocity, the expe- 
rimental heat transfer process can be treated similar to the 
actual heat transfer process in the aggregate hopper. 

Figure 1a shows the sketch of the similarity experi- 
mental system. The cooled air is drawn into the air inlet 
flue of the experimental hopper by the fan, and ex- 
changes heat energy with the stones stored in the hopper, 
finally flows out through the air outlet flue. The aggre- 
gate flow velocity is usually in the order of 0.001m/s in 
actual concrete works, and it is so small that its flow ve- 
locity can be ignored in experiments. So the stones are 
statically accumulated in the experimental hopper. 

The experimental hopper is contracted in equal pro- 
portion according to the actual hopper, and size of it is 1 
mx1 mx2 m (XxYxZ), and size of the air inlet and outlet 
holes is 0.16 mx0.16 m. The actual air average velocity 
in the hopper cross section is about 0.5-1.5 m/s, so the 
converted air volume flow rate is about 3000 m*/h in the 
similarity experiment. As for the fan chosen in the simi- 
larity experimental system, the rated air flow is 2996 
m°/h, the maximum air pressure is 1300Pa, and the rated 
power is 2.2 kW. The flow rate can be adjusted by the 
butterfly valve. Temperature of the cooled air is about 10°C 
lower than that of the stones. The experimental subjects 
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(a) Sketch of the experimental system 
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D (z/L=1) 


C (z/L=2/3) 


B (z/L=1/3) 


A (z/L=0) 


(b) Measuring cross sections 


1. Butterfly valve; 2. Fan; 3. Speedometer; 4. Air inlet flue; 5. Air outlet flue; 6. Experimental aggregate hopper; 
7. Thermocouples; 8. Data acquisition system; 9. Computer. 


Fig. 1 Similarity experimental system 


are marble stones in three sizes: G2 with diameter of 40 
mm-80 mm, G3 with diameter of 20 mm-40 mm, G4 
with diameter of 5 mm-20 mm. The porosity measured 
in experiments is 0.4-0.5. 

Three-dimensional temperature field in the experi- 
mental hopper is measured by K-thermocouples; all data 
are recorded by the data acquisition system and then are 
transferred to a computer. As shown in Figure 1b, ther- 
mocouples are placed on the four cross sections in the 
experimental hopper. Each thermocouple adheres to sur- 
face of one stone by metal foil. There are 12 measure 
points regularly arranged in each cross section. Agilent 
34970A Data Acquisition is used to measure and record 
all temperature data. All the data are simultaneously 
measured and recorded every 10 seconds until the expe- 
riment ends. All thermocouples have been calibrated and 
the absolute error is not more than +1 °C. The pitot tube is 
used to measure the air flow velocity, and its measure 
precision is 0.5%. 


Simulation of the similarity experiment 


The transient heat transfer process in the experiment is 
simulated by means of the porous medium model. There 
are some important assumptions that should be indicated: 

1) Thermal equilibrium between the solid and the gas 
phases is employed in the simulations; 

2) The porous domain is assumed a homogeneous re- 
gion with space-averaged properties; 

3) The gas phase is considered as an incompressible 
ideal gas; 

4) The radiation between particles is ignored. 

The first assumption is suitable for small particles and 
fluids with high heat diffusion [7, 8]. In the aggregate 
hopper, the average air flow velocity in the interstitial 
space is about 1 m/s and the convective heat transfer 
coefficient is about 10 w/m?K [3], the heat diffusion is 


much faster than that of the natural convection or heat 
conduction. In addition, the influence of the stone di- 
ameter can be ignored because the interior heat resistance 
is smaller than the exterior convective heat resistance. 
The Biot number shows the ratio of the interior conduc- 
tive heat resistance and the exterior heat resistance, the 
influence of the diameter can be ignored if Equation (1) 
is valid [9]: 

a:-R/3 4 R (1) 

A 3 

Even for one stone with bigger diameter of 80mm, 
Equation (1) is still valid, so the influence of the diameter 
of the experimental stones can be ignored. 

So, the assumption of the thermal equilibrium adopted 
in the paper is reasonable. A single space-averaged tem- 
perature is then used which can be predicted by a con- 
vection-diffusion equation (2). 


Z lp, +(1-£)p,C, ) +y (p, üC, ) 
=V(k,,VT)+(p,C,) V-(eD" VT) 


where D is an effective diffusion tensor depending on 
fluid velocity and taking into account the axial and radial 
fluid dispersion due to the solid obstacles [10]. 

The transient flow and cooling process in the experi- 
mental hopper is simulated by CFD method, the com- 
bined structure of the stones and the interstitial spaces is 
set as porous domains. About two million tetrahedral 
grids are meshed in the model. The boundary condition 
on the air inlet hole is the flow velocity measured in ex- 
periments; the boundary condition on the air outlet hole 
is outflow. Other input parameters like initial temperature 
of the cooled air, initial temperature of the stones, the 
porosity, are determined according to the experimental data. 

Some typical experimental result and simulation result 
are selected to make the comparison, as Figure 2 shows. 


Bi= 


(2) 
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Fig. 2 Comparison of the temperature field between the experiment and the simulation (After 2100S (G3, 7j=23.4°C, T,=12.8°C, 


uUa=13.8 m/s)) 


The area between two dashed lines in the simulation pic- 
ture corresponds to the measured area. The prediction is 
in good agreement with the experiment. Because of the 
lower air flow resistance in the inlet flue, the cooled air 
easily rushes to the end of the flue and then changes its 
flow direction to disperse into interstitial spaces, so the 
region near the end of the inlet flue are cooler than others. 
Figure 3 compares the measured and predicted tran- 
sient average temperature in two cross sections. Overall, 
the simulations give good results. So, all hypotheses intro- 
duced in CFD simulations are reasonable. The simulation 
method can be reliably used in actual air pre-cooling works. 


Influence sensitivities of working conditions 


Table 1 lists change ranges of working conditions ac- 
cording to actual air precooling works. All 256 permuta- 
tions are simulated by the CFD method mentioned above 
in order to analyze influences of these conditions on 
transient heat transfer characteristics of the aggregate. 
The initial temperature of the aggregate is 293 K, and the 
air inlet temperature is 270 K. 

Fig.4 shows influence of the cross section area of the 
hopper on transient cooling process of the aggregate. The 
larger the cross section area is, the higher the aggregate 
average temperature is. The reason is that the temperature 
field is more nonuniform and the stones far away from 
the air inlet flue are more difficult to be cooled when the 
cross section area becomes larger. 

Fig.5 shows influence of the porosity on transient 
cooling process of the aggregate. When the porosity is 
larger, the stones can more fully touch the cooled air 
flowing in interstitial spaces, so the cooling rate is faster. 
The porosity in the range of 0.4-0.5 has negligible influ- 
ence on the aggregate transient cooling process. 

Fig.6 shows influence of the air flow velocity on tran- 
sient cooling process of the aggregate. The cooling rate 
obviously becomes faster and the average temperature of 
the aggregate significantly decreases when the air flow 
velocity increases. So the air flow velocity has great in- 
fluence on the cooling process. 
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Fig.3 Transient average temperature in cross sections (G3) 


Table 1 Actual working conditions 


Cross section area of 10 18 24 30 
the hopper A, (m7) (2.5x4) (3.6x5) (4x6) (5x6) 


Height between two air 


flues L (m) 3 3.5 4.5 5 
Average air flow velocity 0.2 0.6 1 15 
in cross section u (m/s) 

Porosity € 0.4 0.43 0.46 0.5 


—O— 4,=10 m’, cross section A 


294 
292 
290 
288 


—A— A,=18 m’, cross section A 
—4— A,=24 m’, cross section A 
—O— A,=30 m’, cross section A 


—O— A,=10 m’, cross section C 


2 a —vy— A,=18 m’, cross section C 
b 284 p ? 

E 282 —p>— A24 m’, cross section C 
3 280 —O— A,=30 m’, cross section C 
E 278 

om 


276 
274 
212 
270 


0 2000 4000 6000 8000 10000 
Time (s) 


Fig. 4 Aggregate transient temperature with different cross- 
section area of the hopper (L=5 m, u=1 m/s, €=0.4) 
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Fig.5 Aggregate transient temperature with different porosity 
(L=5 m, u=0.6 m/s, A,=24 m°) 
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Fig. 6 Aggregate transient temperature with different air flow 
velocity (L=3.5 m, A,=18 m?, e=0.4, z/6 =0, cross sec- 
tion A) 


Correlation Analysis of the Transient Heat 
Transfer Characteristics 


The transient heat transfer regularity of the aggregate 
can be deduced according to the enormous simulation 
results under different working conditions. Nomogram is 
an available and convenient method to analyze the tran- 
sient heat transfer process and is commonly used in en- 
gineering. Nomogram of the air precooling process of the 
aggregate is drawn in the paper, and the fitting correla- 
tions are suggested. 


Equivalent thermophysical properties of the porous 
structure 


The structure composed by the stones and the intersti- 
tial spaces in the aggregate hopper is treated as a porous 
structure. The equivalent thermal conductivity 4, is com- 
monly determined by fluid and solid conductivity, poros- 
ity and solid structure, while the fluid convection in in- 
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terstitial spaces is usually ignored [5, 11, 12]. However, 
the air convection between stones is so critical in the air 
precooling process that it couldn’t be ignored. 

Fig.7 shows assumed regular arrangement of two ad- 
jacent stones. It is assumed that the heat energy transfers 
from the lower surface to the upper surface of the unit. 
The heat transfer capacity in the unit includes two parts: 
one is the convection heat transfer in the interstitial space; 
the other is the heat conduction between two adjacent 
stones. 


Fig.7 Arrangement unit 


The convection heat transfer in the interstitial space 
can be calculated by Eq. (3): 


ee AT -2axdx _mdA,-AT - (20d 
2.) ( d/2Y =a 
K (3) 
mA? -AT| | d ar) A A 
+—* | In| —,/1 -| — | +— |-In| —— 
An. 2 d An Any 


Where AT is temperature difference between the lower 
and upper surfaces of the unit. 
The heat conduction between two adjacent stones can 
be calculated by Eq. a 
nr -AT  ar*Ag- AT 
Q,= = — (4) 
H HAs Ay H 
So the equivalent thermal conductivity of 2, can be 
deduced from Eq. (3) and Eq. (4): 


2 
= A+ 274 hora) +2 
aAT/H 2d apd 


d ary A A, \|. ara 
In} —,j1-—] — | +— |-ln| — | |+ z 
2 d) aM Op, d? 
The diameter of the touch area has little change and 


always closes to zero in the porosity range of 0.4-0.5. So 
Eq. (5) can be simplified as follows: 


2 
gs In d A —Iln As (6) 
2 ad 2 Oy Ap. 


e 
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It is assumed that the air flow in each interstitial space 
is one-way from lower cross section to upper cross sec- 
tion in the hopper, so the equivalent diameter of the as- 
sumed flow pipe can be calculated by Eq. (7): 


A(d?-nd?/4) (4 
de (F} 4 4 


The convective heat transfer coefficient in the as- 
sumed flow pipe is given by: 


A 
Ot, =0.023:—“ . Re?® Pr °* (8) 
kx 
Re, = Lalli (9) 
Ha 
u = = (10) 
E 


The equivalent thermal conductivity of 4e can be cal- 
culated by Eq. (6) — Eq. (10). 

The equivalent density and specific heat of the com- 
bined structure are determined by Eqs. (11) - (12): 


P eP 1-2), (11) 
a lna g E (12) 
Pe Pe 


Nomograms abstracted from simulation results 


In air precooling aggregate works, four side walls of 
the aggregate hopper are adiabatic and heat can only be 
transferred from lower surface to upper surface in z di- 
rection. So the transient heat transfer process of the ag- 
gregate can be approximately described by Eq. (13) [9]: 


o(zt)_T(z,t)-T, (Fo Bi, z) (13) 


bo To Z Tı 
3 aT 
The Fourier number: Fo= F (14) 
, _ ao 
The Biot number: Bi= ra (15) 
P A 
The heat diffusivity: a=— (16) 
pc 


The cooled air flows in the inlet flue and then dis- 
perses into interstitial spaces in the hopper, which is dif- 
ferent from the ideal convective condition on the upper 
surface and the lower surface of the hopper. So it is ne- 
cessary to redefine the Biot number and the Fourier 
number. The effective Biot number is redefined accord- 
ing to the converted convective condition in the whole 
cross section (see Eq. (17)); the effective Fourier number 
is redefined according to the hopper size and thermal 
properties of the combined structure (see Eq. (18)): 


OE 


Bi (17) 


Fo=— (18) 
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Other parameters are defined as follows: 
The effective heat diffusivity: 


a; = Ae. (19) 
P e Ce 
The virtual convective coefficient: 
a, =0.023- Me Re? Pr °* (20) 
D, 
Equivalent diameter of the hopper: 
4A 
D, = U (2 1) 
Virtual Reynolds number: 
D 
Re, = Pak Pe (22) 
La 


So the transient heat transfer process in the aggregate 
hopper can be described by Eq. (13), Eq. (17) — Eq. (22). 
Calculation of Eq. (13) is so difficult that nomogram is 
usually introduced to describe and predict the transient 
heat transfer process, and it has been widely used in en- 
gineering. X-coordinate of the nomogram is Fo; Y-coor- 


-1 


dinate is . In the paper, the transient heat transfer 


a 0 
characteristics of the aggregate are analyzed according to 
the enormous simulation results, Bi and Fo are calculated 
by Eqs. (17) — Eq. (22), and then nomograms under 
different work conditions are drawn. 

Simulation results show that the transient heat transfer 
rule is similar when L/D, changes little. So the simula- 
tions are ranked according to L/D, (see Table 2). 

Nomograms are selectively shown in the paper. Fig. 8 
show nomograms in cross section D (z/L=1). 


Table 2 Ranked simulations according to L/D, 


Cross section Height 


No. L/D, LID. area of the between two 
hopper 4, (mô air flues L (m) 

0.625 24 3 
1 0.634 

0.642 30 3.5 

0.717 18 3 
2 0.724 

0.73 24 3.5 

0.826 30 4.5 
3 0.831 

0.836 18 3.5 

0.917 30 5 
4 0.943 0.938 24 4.5 

0.975 10 3 

1.04 24 5 
5 1.058 

1.075 18 4.5 

1.137 10 3.5 
6 1.166 

1.194 18 5 

1.46 10 4.5 
7 1.543 

1.625 10 5 
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Fig.8 Nomograms when L/D, =1.543 (cross section D) 


Nomograms show that the compared excess tempera- 


-h 


ture increases when Bi increases, which means 


a 0 
that temperature of the aggregate more obviously de- 
viates from the initial value and approaches air tempera- 
ture when the air convection becomes intensive. 


Fitting equations for transient heat transfer processes 


-h 


Change regularity of seriously depends on 


a 0 
Fo and Bi, as shown in nomograms. The exponential 
function of Fo can be used to predict the compared 
excess temperature according to the obtained nomo- 
grams: 
T-T 
T- 


Where A, B, C, Fo* are fitting coefficients and se- 
riously depends on Bi. Fo* explains temperature delay of 
the aggregate in the cross section. In the cross section A 
(z/L=0), the aggregate temperature immediately changes 
when the air flows into the inlet flue, so Fo*=0. 


= A-B-exp|-1-C-(Fo-Fo*)] (23) 


Table 3 Correlations of the fitting coefficients 


Cross section 


Correlations of the fitting coefficients can be deduced 
from simulation results, as shown in Eq. (24) and Table 3. 
A=A" BitA"=f (L/D.,Bi) 
B=B' -Bi+B"=f (L/D, Bi) 
C=C" -Bi+C"=f (L/D, Bi) 
*=f (Bi) 

Fig. 9 presents relative error of the fitting equations 
for the compared excess temperature in four cross sec- 
tions. The relative error is about less than +10%, so the 
fitting equations are precise enough to predict the tran- 
sient heat transfer process of the aggregate in the air pre- 
cooling hopper. 


(24) 


Conclusions 


The transient heat transfer characteristics of the ag- 
gregate in the air precooling hopper are analyzed by CFD 
method, nomograms of the transient heat transfer 
processes are sorted out and fitting equations are deduced 
to precisely predict the transient heat transfer process. 
Some important conclusions can be obtained from the 
analyses: 


Fitting coefficients 


C'=1366.4- L/D, + 554.02 ; C" = -530.2 - L/D, + 215.72 ; Fo*=0.0005: Bi-1.166 


C'=1126.9- L/D, + 489.25 ; C" = —621.13- L/D, + 253.8 ; Fo*=0.001:Bi-1.146 


A (z/L=0) A' = -0.0225 L/D, + 0.0453 ; A" =0.0451- L/D, +0.8156 ; 
B' =-0.0254- L/D, +0.0192 ; B" =0.0919- L/D, +0.5492 ; 
C' =1332.8- L/D, +154.81 ; C" =-275.36- L/D, +114.87 ; Fo*=0 
B (2/L=1/3) A'=-0.0123- L/D, +0.0219 ; A" = 0.0535- L/D, + 0.807 ; 
B' =0.0037- L/D, —0.0351 ; B" =0.5833- L/D, +0.4248 ; 
C (z/L=2/3) A' =—0.0838- L/D, + 0.0246 ; A" =0.2522- L/D, +0.7434 ; 
B' =0.0374- L/D, —0.1893 ; B" =1.0151- L/D, + 0.2248 ; 
D (z/L=1) A' = —0.0838- L/D, +0.0246 ; A" =0.2522- L/D, + 0.7434 ; 


B' =0.0374- L/D, —0.1893; B"=1.0151- L/D, +0.2248 


C'=1126.9- L/D, + 489.25; C"= 


—€21.13- L/D, + 253.8 ; Fo*=0.0016- Bi-1.14 
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Calculation of fitting equations 
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Simulation result 


(a) Cross section A (z/L=0) 


A A,=10 m? 
o A,=18 m? 
O A,=24 m? 
y A,=30 m? 


Calculation of fitting equation 


0.0 0.2 0.4 0.6 0.8 1.0 
Simulation result 
(c) Cross section C (z/L=2/3) 


Calculation of fitting equation 


Calculation of fitting equation 
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A A,=10 m? 
o A=18 m? 


O A24 m? 


0.0 0.2 0.4 0.6 0.8 1.0 
Simulation result 


(b) Cross section B (z/L=1/3) 


0.0 0.2 0.4 0.6 0.8 
Simulation result 
(d) Cross section D (z/L=1) 


Fig.9 Relative error of the fitting equations 


1. The structure composed by the aggregate and the 
interstitial space in the hopper can be treated as a porous 
structure, and the similarity experiment results verify that 
the space-average method is suitable for the stones with 
diameter not more than 80 mm. 

2. Size of the hopper and the average air velocity in 
the cross section have obvious influence on the transient 
heat transfer process of the aggregate, while the porosity 
in the range of 0.4-0.5 has little influence, according, to 
enormous simulations of actual air precooling works. 

3. Correlations of the compared excess temperature 


-h .. , . 
7 fitted from nomograms can precisely predict 
a *0 

the air precooling transient heat transfer process of the 

aggregate. The calculation offers a reliable method for 

design of the air precooling process to avoid temperature 


cracks and unreasonable energy consumption. 
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